The simplest model of polymer crystal exhibiting polymorphism. 
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Almost all the polymer crystals have several polymorphic modifications. Their structure and 
existence conditions, as well as transitions between them are not understood even in the case of 
the 'model' polymer polyethylene (PE). For analysis of polymorphism in polymer crystals, we con- 
sider the simplest possible model of polymer chain: an extended fiat zigzag made of 'united' atoms 
(replacing ClU-groups in PE chain); the united atoms belonging to different zigzags interact via 
Lennard- Jones potential. Analysis of potential of interaction between such zigzags allowed to pre- 
dict the structure of five possible equilibrium lattices in polymer crystal built out of such zigzags. 
Molecular dynamics simulation of the crystal built out of flexible zigzags showed that, depending on 
model parameters (dimensions of the zigzag and equilibrium distance of Lennard- Jones potential), 
one to three of these lattices are stable in bulk at low temperatures. We have determined the model 
parameters at which the existing stable lattices are analogous to the ones observed in real PE and 
linear alkanes. The triclinic lattice has the lowest potential energy, then follows the monoclinic lat- 
tice, and the orthorhombic lattice has the highest energy, exactly as in full atomic (with hydrogens) 
molecular dynamics models. 



INTRODUCTION 

Polymorphism is an attribute of polymer crystals. It 
is pronounced even in the simplest, 'model' polymer 
polyethylene (PE) [-Cfi2-]„. The stable low-temperature 
phase in its linear odd n-alkanes is the orthorhombic (O) 
one, in its even n-alkancs - the triclinic (T) one up to 
n=24 and the monoclinic (M) one when 26 < n < 34 [1] . 
PE melt and solution crystallization gives O folded-chain 
crystallites, while on substrates [2, 3] and after polymer- 
ization inside nanochannels [4] there is a large amount of 
non-0 (M and T) crystallites. At higher temperatures, 
one can also see a series of rotator phases in n-alkanes [5] , 
and the hexagonal conformationally disordered phase in 
PE [6]. 

There exists the long history of study of PE phases de- 
pending on pressure and temperature. Starting with the 
works by Seto [7], the transitions between orthorhombic 
and monoclinic phases were widely investigated (see, for 
example, [8-10] and references in them; for a review, see 
[11]). On the other hand, the P-T phase diagram near the 
triple point melt-orthorhombic-hexagonal phases of PE 
was very early obtained [12] and then many times speci- 
fied depending on length of PE chains, samples morphol- 
ogy and thermal history. But, in spite of all the efforts, 
the exact place of M and T phases on this diagram has 
never been established; the character of their instability 
is still unclear. Almost at any pressure, they normally 
coexist with the O phase. The pressure, at which they 
appear and the rate at which they transform back to the 
O phase, depend on mode of preparation of the samples 
(compare, for example, [7] and [13]). The understanding 
of the structure of these lattices can shed light on this 
unsolved problem. 

One usually treats polymorphic modifications in 



molecular crystals in the framework of the 'closest pack- 
ing' approach [14], considering number of 'contacts' be- 
tween atoms in different lattices. This approach, as well 
as molecular mechanics studies [15], does not allow to 
say which interactions between molecules in a lattice are 
principal and what balance of forces forms every lattice. 
Correspondingly, one can not say which influence on the 
lattice is able to destroy its system of couplings and to 
provide the transition to another lattice. For example, 
in PE, T and M phases appear under pressure [7, 13], 
tension [16] and shear deformation[17]. But the needed 
directions of influence and resulting modes of deforma- 
tion are still not understood. 

To analyze a coupling between two molecules in a 
molecular crystal, one is to have their potential of inter- 
action in treatable (better analytical) form. In the case 
of PE, it is natural to suppose that the details of the 
potential determined by hydrogen atoms are not princi- 
pal for the system of couplings in the lattice. Then one 
can replace CH 2 -groups by 'united' atoms (UA), and the 
chains in the lattice - by flat extended zigzags made of 
the UA. 

However, it is a common opinion probably since 1980th 
[18] that 'the UA approximation docs not provide a good 
representation of the crystal structures of alkanes' [19]. 
Indeed, we saw [20, 21] that, at some model parameters, 
the only crystal phase at low temperatures is a M lat- 
tice, different from the M lattice of PE. The other model 
parameters can even give a mobile hexagonal phase at 
300K [19]. Therefore, all the works on crystalline phases 
in PE and alkanes use the full atomic (FA) model (with 
hydrogens) in spite of 'an order of magnitude increase in 
the processor time compared to the UA case' [19]. 

The goal of this article is to analyze equilibrium crys- 
tal lattices present in UA PE model at different model 



2 



parameters, compare them with the lattices present in a 
realistic FA model of PE (taking into account hydrogen 
atoms) and to find the parameters at which the stable 
in bulk at low temperatures lattices are analogous to the 
ones observed in the FA model - and in real PE and linear 
alkanes. 

The FA model that we will use as a standard was intro- 
duced in [22]. It uses generally accepted Amber99 force 
field [23] and gives the results on energies and densities of 
PE lattices similar to the ones obtained in other realistic 
empirical force fields [15, 24]. 



GEOMETRIC PARAMETERS OF UA PE 
MODEL. 

UA PE model is widely used in modelling of properties 
of liquid and gaseous n-alkanes. The review [25] com- 
pares different sets of constants for simulation of fluid 
properties. We shall dwell on geometric parameters re- 
sponsible for the presence of polymorphism in crystalline 
phase. 

Van-der-Waals interaction of two UA belonging to dif- 
ferent zigzags and placed r apart is normally described 
by 6-12 Lcnnard- Jones potential: 
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(1) 



The first evident geometric parameter of the model is 
the degree of 'collectivity' of interactions between UA, 
the ratio ro = Rq/c, where c is longitudinal period of 
the PE chain. The more is ro, the more UA of the 
neighboring zigzags take part in generation of the po- 
tential near an atom [26]. For example, in the Ambcr99 
force field [23], for hydrogen- hydrogen, hydrogen-carbon 
and carbon-carbon interactions, the values r equal to 
1.166, 1.331 and 1.497 correspondingly. This parameter 
is widely varied in different simulations. 

The second geometric parameter is the same in all but 
one of the works. Namely, UA are as a rule placed on car- 
bon atoms, although, intuitively, it seems more adequate 
to shift united force centers closer to the hydrogens in the 
CH 2 -groups. It has been done in model [27] where mass 
and force centers were divided (mass centers stayed on 
carbon atoms). According to work [25], this model gives 
the best agreement with experimental data in modelling 
of fluid properties of n- alkanes. We did not divide force 
and mass centers of UA, and placed a united atom on the 
bisector of the angle between hydrogen atoms, keeping 
the longitudinal period c of the PE chain and changing 
the distance d between the two rows of UA (fig. 1). In 
such an approach, we explicitly introduce the natural ge- 
ometric characteristic of a flat zigzag, the ratio do = d/c, 
as the second geometric parameter of the model. When 
a united atom shifts from carbon to the point between 
hydrogens, do changes from 0.337 to 0.843. 




FIG. 1: Choice of the UA' positions relative to the PE chain. 
Carbon (large) and hydrogen (small) atoms are white, UA are 
black, the bonds between the UA are dashed. 



EXTREMA OF POTENTIAL OF INTERACTION 
BETWEEN TWO EXTENDED PARALLEL 
ZIGZAGS OF UA. 

If we accept that two UA interact according to (1), 
then, if ro is reasonably large, ro > 1.115, the potential 
of interaction of one united atom with an infinite row of 
UA can be obtained by Poisson summation formula (see 
fig. 2(1)) where the first two terms quite suffice [26]: 



Jo{ko,(i) + Ji(fc ,/3)cos (C) + 



Wi(*0>AO = j^k e +J 2 (fc ,/3)cos(2C) + ...+ 
+J m (k ,P)cos{mC) + ... 

where ko = 27r?*o, /? = 27r6/c, £ = 2ttz/c, 
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The function Jo(ko,P) has a minimum equal to 
-1.43132... at the point /3 = O = B k , B = 
0.9471316..., and the function J\ (ko, P) can be very good 
approximated by an exponent near the point f3o'- 

J\ (ko,P) « A x (k ) cxp (-A(fc ) (0 - A,)) , 



Ax (k ) = -k 5 exp(-Po) fx(Po,k ) , 



A(fco) = 1 - 



dfi(P,k ) 



dp 



fi (j3o,k ), 
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FIG. 2: Steps of analysis of potential of interaction between 
two infinite parallel extended zigzags made of UA. (1) Find- 
ing potential Wi of an infinite row of UA at the point P as a 
function of b and z. (2) Finding potential of interaction Wa 
between two zigzags. We show the projection of the zigzags 
on a plane orthogonal to their axes. The arrow points from 
the row of even atoms to the row of odd atoms of one zigzag. 
Black circle denotes the row in which zeroth atom belongs to 
the projection plane, white circle - the row in which the near- 
est to the projection plane atoms are at the distance c/2 from 
it. We find W$ as a function of the distance D between the 
zigzags, their setting angles 9\ and 02, and the shift z of the 
second zigzag along its axis from the projection plane. When 
z = 0, the rows of even and odd atoms of the second zigzag 
also become 'black' and 'white' correspondingly. (3) Extrema 
of W4. In the sketches, we show only the arrows introduced 
on step (2); here, we imply that the rows of atoms are 'black' 
and 'white'. The lowest in energy coupling (the absolute min- 
imum) has notation A, the higher extrema - notation B, C,... 
and so on. 

fi(P,ko) = 

6 i / 945 + 945/3 + 420/3 2 + \ _ 
- fe o iso^TT ^ +10 5/3 3 + 15^ 4 + /3 5 ) 
-2 ((3 2 +3/3 + 3) //3 5 

To obtain the energy of interaction Wi between two 
zigzags (per 1 united atom), one is (see fig. 2(2)) to sum 
four terms of the type W\ (2), namely, the energy of 
two successive UA of one zigzag in the presence of two 
rows of another zigzag (and to divide by four). Unfor- 
tunately, only in the case of proximity of the zigzags to 
round cylinders d/D ~ do/r « 1 the function W4 can 
be adequately presented by the first terms of its Fourier 
series in angles 6*i and 8 2 . This condition does not hold 
in models exhibiting polymorphism in crystal lattices. 
Therefore, we found the extrema numerically as a sum of 
the mentioned four terms. 

It is evident that £ = or 7T for equilibrium config- 
urations. Taking £ = for dcfinitcness sake, we ob- 
tained four (physically different) stationary points (see 
fig. 2(3)): A (0r = tt/2, 9 2 = -tt/2), B (6 1 = tt/2, 
6 2 = tt/2), D (0i = 0, 6 2 = 0), E (6»i = 0, 9 2 = tt) similar 
to the ones in dipole-dipole interaction between electrical 
or magnetic dipolcs. There are four more couplings des- 
ignated as C in fig. 2(3). We did not draw arrows here, 
because all four combinations are possible, although the 
setting angles are a little different. 



TABLE I: Distances ys = Ds/Ro and energies ws = W4/E 
of equilibrium configurations of two zigzags (extrema of W4) 
listed in fig. 2(3)). The presented numerical values are calcu- 
lated for model parameters G4 (do — 0.647, ro = 1. 541). In 
other cases, the extrema and their order are the same, the 
numerical values slightly differ. 
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1.117 


-2.599 
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-3.681 
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0.912 


-3.894 



Among the listed couplings, A is a stable minimum, 
the rest stationary points are saddles. B is unstable in z 
direction (slides into A), C - in 9 2 direction (slides into 
B or A), D - in any direction in (9\, 9 2 ) plane, as well 
as E, which is also unstable in z direction. The equilib- 
rium distances yg — Dg/R and energies w$ = W4/E 
(E(ro,e) = 1.6862roe) in these configurations are listed 
in Table I. The unit of energy E (ro, e) is chosen so that 
the first term in (2) has minimum equal to (-1). 

FORMATION OF EQUILIBRIUM LATTICES. 

We accept that structure of lattices is basically deter- 
mined by interactions of a zigzag with 6 (or may be 8) 
zigzags of its first coordination sphere. Let us sketch a 
lattice as a set of arrows (introduced in fig. 2) represent- 
ing zigzags in the projection plane orthogonal to the axes 
of the zigzags (see fig. 3). 

Lattices based on A and B couplings: T and Mi. 

Geometrically, it is impossible to build a lattice using 
only the strongest and stable A couplings. One can form 
a row of arrows (for example, horizontal) so that every 
two neighboring arrows point in opposite directions (one 
- up, the second - down), and this will be a stable struc- 
ture; but one still is to build a two-dimensional pattern 
from these rows. Combining the rows, one can not avoid 
forming couplings of D and E types between molecules 
in neighboring rows. As a result, we get a T lattice. 

In it, A couplings in horizontal rows store about twice 
energy than the couplings between the rows (compare the 
energies of A, D and E couplings in Table I) and prevent 
the zigzags from rotation. E couplings are unstable re- 
garding the shift of the zigzags along their axes, and this 
instability is compensated by the others couplings. 

T lattice built out of rows with A couplings and having 
couplings of D and E types between rows may happen to 
be energetically less favorable than lattice Mi built out 
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FIG. 3: Sketches of possible types of equilibrium lattices built 
out of zigzags interacting through potential having extrema 
presented in fig. 2 (3). Every lattice consists of (horizontal) 
rows of zigzags connected by equilibrium couplings: T - A, 
Mi - B, O - C, M2 - D. In real phases (see fig. 4), not only 
the couplings between rows but also the couplings forming 
rows may be slightly out of equilibria. 



of rows with B couplings but having only D couplings 
between the rows. 

A separate row with B couplings is an equilibrium 
structure, but it is unstable regarding the shift of any 
zigzag along its axis. In lattice M 1; the rows with B cou- 
plings stabilize by cohesion with each other: the arising 
couplings of D type are stable in z direction and com- 
pensate the instability of B couplings. 

In both lattices T and Mi, the couplings D and E be- 
tween rows are off their equilibria (although the forces 
are small), and every zigzag needs two couplings with 
the neighboring row to keep its setting angle. 



Lattices using C and D couplings: O and M2. 

In lattice M2, horizontal D couplings in rows do not 
form this lattice but rather are its vulnerability. The 
strongest couplings are the non-equilibrium ones between 
rows, intermediate between B and D, and designated as 
BD in fig. 4. They appear in both lattices M2 and O. 
Again, to maintain the equilibrium, every zigzag needs 
two different BD couplings with another row. In lattice 
M2 , the setting angles of zigzags in D couplings may be 
off their equilibria, and BD couplings 0-1 and 1-2 may 
have different lengths, and so one of these couplings may 
be closer to B than another. 

There can be at least two different equilibrium O lat- 
tices. In both of them, the couplings 0-2 and 1-2 (fig. 4, 
at the bottom) are equivalent and therefore the sum of 
the setting angles of zigzags Q\ and 62 is tt. In lattice Oi, 
6» 2 - 6>i < tt/2, in lattice 2 , 2 - 0i > 7t/2. 




FIG. 4: Lattices T, M2 and O2 (viewed in projection on a 
plane orthogonal to chains' axes) in UA PE model (G4: do = 
0.647, ro = 1.541) in comparison with analogous phases of 
crystalline PE (FA MD simulation, the details of the model 
are in [22]). Radii of the balls are equal to van-der-Waals radii 
of the corresponding atoms (carbon, hydrogen and 'united' 
ones). 



As we remember, the coupling C is unstable regard- 
ing rotation of zigzag 2 (see fig. 2) and stable regarding 
rotation of zigzag 1. In accord with it, in O lattice, the 
zigzag is prevented from rotation basically by its neigh- 
bors for which arrows are aligned along the direction to 
the zigzag. Two other C couplings give almost flat po- 
tential, and BD couplings only shift a little the minimum 
of the energy. 



CHOICE OF MODEL PARAMETERS FOR PE. 

To find the point the best corresponding to PE in the 
rectangular area of reasonable parameters, we selected 
140 points in the rectangle (intersections of grid lines in 
fig. 5) and studied possible equilibrium lattices in every 
point. 

In the computational cell, we put 4 rigid zigzags par- 
allel to z axis (more exactly, two UA from every of 4 
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FIG. 5: Domains of existence of different lattices depending 
on model parameters (for details, see section ). The points 
the best corresponding to real PE lie in the white area of the 
domain bounded by dashed polygon; the point G4 is our best 
choice. 



zigzags). Every such zigzag has 4 coordinates, 3 spacial 
ones and the setting angle as introduced in fig. 2. The 
first zigzag was placed so that its spacial coordinates were 
zeroes, and so we got 13 variables. Two edges of the com- 
putational cell were parallel to x and z axes, and so we 
got 3 variables characterizing the computational cell (two 
lengths of the edges in (x, y) plane and the angle between 
them). The computational cell was repeated nine times 
in directions of its edges in (x, y) plane, and 80 times in 
z direction. In all three directions, we imposed periodic 
boundary conditions. For the full Lennard- Jones energy 
of the obtained crystal, we solved the minimization prob- 
lem in the space of 16 variables starting at structures Mi, 
M2, Oi, O2 and T. In points G3 and G2, we used the lat- 
tice parameters obtained in approximation of the first 
coordination sphere with the use of analytical potential 
(2). After we got the solutions in these points, we moved 
to the next points using at every step as initial structures 
the ones obtained in the previous points. The results are 
presented in fig. 5. 

T lattice exists almost everywhere in the rectangle, ex- 
cluding a couple of lowest rows of points. In the 10% grey 
domain 'distorted T', lattice T is distorted, the setting 
angles of the zigzags differ from 7r/2 and — 7r/2, contrary 
to the experiment [7], and so possible candidates for PE 
can not lie in the domain 'distorted T'. 



TABLE II: Lengths Y>tj of the couplings (in angstroms) and 
setting angles of the zigzags (in degrees) in lattices T, M2 and 
O2 (fig. 4) at different model parameters (in points shown in 
fig. 5) in comparison with FA model. In the brackets, we 
point out the type of the coupling. 
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Lattice O2, corresponding to O lattice of real PE, exists 
everywhere except for 60% grey domain designated as 
'no O2'. The setting angle Q\ changes from ~ 30° in the 
upper part of the rectangle to 45° at the boundary with 
'no O2' domain. Lattice Oi is present everywhere, and 
has lower energy than lattice O2. 

Lattice M2, corresponding to M lattice of real PE, ex- 
ists everywhere except for the 40% grey domain desig- 
nated as 'Mi' where the only minimum of M type is lat- 
tice Mi. In the hatched domain 'M2+M1' these lattices 
coexist, and lattice M2 has lower energy. 

So, in the rectangle, only the points not painted in grey 
can be possible candidates for modelling of PE. Among 
them, the points belonging to the line G^-Gq have the 
density equal to the density of PE at low temperatures. 
So, it looks reasonable to select the point for PE in the 
domain bounded by a dashed polygon. 

To make more precise choice, we compared the sizes 
of the lattices in several points on the line G^-Gq and in 
the FA MD model. The results arc presented in Table 
II. One can see that the proportions of the lattices in 
toto approach the ones in FA model when we move from 
the point Gq to the point Gt- Besides, the dynamical 
behavior of the crystal will probably be closer to the real 
one in the points with less d because mass centers of UA 
will be closer to carbons as it is in real PE. 

The chosen dashed polygon in fig. 5 is divided into 
two principally different zones. On the left, the models 
have four different lattices: T (ground state), M2, Oi 
and O2; and Oi was never found in real PE or alka- 
nes. On the right, there is even five lattices: T (ground 
state again), M2, Oi, Mi and O2, where two 'superflu- 
ous' lattices Oi and Mi are again lower in energy than 
the widespread at high temperatures lattice O2. The 
natural question is: are the found lattices really present 
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as phases in the model or are they somehow unstable? 
Indeed, our minimization problem was solved only in 16- 
dimensional space, and we definitely could block some 
modes of transformation between lattices, and some of 
the lattices may prove to be saddles, and not minima. 

STABILITY OF EQUILIBRIUM LATTICES IN 
UA PE MODEL - MD SIMULATION. 

The quickest way to answer the question put in the 
previous section is to carry out a MD simulation of all 
these lattices in large enough sample, in the framework of 
UA PE model having corresponding geometrical param- 
eters, but with normal flexible zigzags, having all mobile 
UA. We took as a computational cell a parallelepiped 
having 120 zigzags 80 UA long packed into expected lat- 
tices. Because we are interested in stability in bulk, we 
imposed periodic boundary conditions in all three direc- 
tions (and so the zigzags became 'infinite'). But, because 
we are interested in the easiest conversions between lat- 
tices, we opened a way for transformations closed in real 
experiments. Namely, we allowed the computational cell 
to change its shape (3 lengths and 3 angles) through uni- 
form change of distances between atoms - according to 
the arising stresses (local stresses averaged over the sam- 
ple). So, wc facilitated conversions between lattices close 
in sizes - if the sample is able to preserve momentum 
and angular momentum during the conversion. In real 
experiments, martensitic transitions in bulk normally go 
through splitting samples into domains because this way 
of transformation is closed. 

The exact values of force constants one uses in the 
model (the value e in Lennard- Jones potential (1), rigidi- 
ties of bonds between UA, of valence and torsional angles 
in zigzags) are not really important when one studies be- 
havior of lattices at low temperatures. We carried out the 
MD experiments at T=1K, and used force constants close 
to the ones chosen by Toxvaerd [27]. Our force constants 
are listed in [21], except for the value e, which some- 
times was taken to be equal to 0.16 kcal/mol instead of 
0.12 kcal/mol. To maintain the temperature, a collisional 
thermostat [28, 29] was used with the parameters mo = 
1 a.m.u. and A = 5.5ps _1 , and a Berendsen barostat 
[30] was used for maintaining pressure in 1 atmosphere. 

We carried out MD simulations of all obtained in the 
previous section lattices at model parameters correspond- 
ing to typical points (having different sets of equilibrium 
lattices) in fig. 5: Gi, G2, G3, Gt and G4. 

One gets the point Gi (d = 0.337, r = 1.673) if one 
puts UA on carbon atoms and then chooses i?o so that, at 
temperature T=300K, the density of the resultant crystal 
(Mi) coincides with the density of O lattice of PE. In 
Gi, we expected two lattices, Mi (ground state) and Oi 
(which is present everywhere). Lattice Mi appeared to be 
stable, Oi transformed into ground state Mi. The mode 



of transformation is very simple and evident; the setting 
angle 6\ (see fig. 4) grows up to 90° and the setting angle 
02 drops down to 90°, the sample slightly changes its sizes 
and we get an accurate Mi lattice. This transformation 
can so easily take place because the rotation of a couple of 
zigzags exactly conserves angular momentum. So, in the 
point Gi, there is only one stable phase, as was already 
known [20]. 

In the point G 2 (d Q = 0.337, r = 1.273), with UA 
placed on carbons but having more 'local' interactions 
between atoms than in the point G\, we expected the 
lattices M2 and again 0\. Both of them proved to be 
stable, and so, in the point G2, we have two stable phases. 

In the point G3 (d = 0.493, r = 1.528), we expected 
Mi, T, Oi and O2. Oi transforms (exactly as in the 
point Gi) into the ground state Mi, the rest lattices are 
stable. So, in point the G3, we have 3 stable phases. 

In the point Gt (do = 0.584, r = 1.549), correspond- 
ing to the choice made by Toxvaerd [27] for liquid alkanes, 
we expected all the possible lattices in the order T, M2, 
Oi, Mi and 2 . Only two lattices, the lowest (T) and 
the highest (O2) in energy proved to be stable phases, the 
rest lattices transformed into the ground state T. We can 
conclude that the domain 'M2+M1' should be excluded 
from the area of possible candidates for PE. 

In the point G4 {do = 0.647, ro = 1.541), we expected 
T, M2, Oi and 02- Only lattice Oi transformed into T; 
the rest lattices are stable. 

So, in the point G4, we have three low-temperature 
phases exactly analogous to the phases observed in real 
PE (fig. 4). The ground state of the system is T phase, 
followed by M2 and then O2, exactly as in the FA PE 
model. The densities of the phases differ from the ones 
in the FA PE model less than by 1%. 

CONCLUSION. 

Potential of interaction of two flat parallel rigid (in- 
finite) zigzags built out of UA has 5 different extrcma. 
Connecting zigzags by couplings close to these extrema, 
one can build 5 different equilibrium crystal lattices, in- 
cluding lattices T, M2 and 2 observed in PE and alkanes 
at low temperatures. So, the presence of hydrogens is not 
necessary for the presence of polymorphism in this poly- 
mer crystal; UA PE model catches the main features of 
interaction between zigzags and provides all the needed 
polymorphic modifications in the crystal built out of such 
zigzags. However, these equilibrium lattices can be sad- 
dles as well as local minima in potential energy of the 
crystal lattice depending on geometrical parameters of 
the model. In different domains of parameter space, one 
to five of these lattices can be local minima. 

Unfortunately, the usual choice of UA positions on car- 
bons eliminates T lattice as a minimum and, at (also usu- 
ally used) very long-range interactions between UA, also 



7 



eliminates O lattice. This is the reason why UA model 
has been in disrepute. But the shift of UA to hydrogens 
leads to appearance of minima for both lattices T and 

o 2 . 

However, a lattice can be a minimum in the space of 
variables of possible crystal lattices having corresponding 
symmetry restrictions, but to be a saddle in full space of 
variables, because the symmetry restrictions may block 
the modes of transformation between lattices. The sim- 
plest way to check if an equilibrium lattice is a phase 
(meaning stable in bulk at low temperatures) is to carry 
out a MD simulation of it in a large enough sample, with 
all mobile UA. We did it for several typical points in space 
of model parameters and found that UA model can have 
one to three phases. These three phases may be (the 
point G4) or not be (the point G3) analogous to phases 
in PE. This narrowed the area for possible candidates for 
modelling of PE to the white part of the domain bounded 
by dashed polygon in fig. 5. Based on analysis carried 
out up to now, we would say that the best choice for 
PE is point G4, although, after the further investigation 
of temperature behavior and phase transitions, the other 
points in the area can prove to be more suitable. 

In the framework of UA model, one can observe phase 
transitions in polymer crystals in reasonable time of cal- 
culations, which makes this model (having properly cho- 
sen parameters) a useful instrument for studying features 
of polymorphic modifications and martensitic transitions 
in polymer crystals. 
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